Hybrid Goldstone modes in multiferroics 
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By using polarized inelastic neutron scattering measurements, we show that the spin-lattice quan- 
tum entanglement in mutliferroics results in hybrid elementary excitations, involving spin and lattice 
degrees of freedom. These excitations can be considered as multiferroic Godstone modes. We argue 
that the Dzyaloshinskii-Moriya interaction could be at the origin of this hybridization. 



Spontaneous symmetry breaking is a powerful concept 
at the basis of many developments in physics, especially 
in condensed matter and high energy physics. The low 
symmetry phase is described by an order parameter asso- 
ciated with low energy and long wavelength excitations, 
restoring the original high temperature phase symme- 
try. These Goldstone modes [l[ are nothing but longi- 
tudinal and transverse phonons in solids, or spin waves 
in magnets. In case of multiferroic materials, two or- 
der parameters, namely the ferroelectric polarization and 
the magnetization, coexist and are strongly coupled by a 
spin-lattice interaction As a result of this entan- 

glement, the multiferroic Goldstone modes are expected 
to be new spin and lattice hybrid excitations called elec- 
tromagnons 0, 0, S 01 ■ While their existence has been 
theoretically predicted for a long time 8, 9], their dual 
nature as both spin and lattice excitations makes them 
challenging to observe and study experimentally. Due 
to their dipole electric activity, they first could be de- 
tected by optical measurements [1, [l, |^ : evidence for 
their existence has been recently reported in different 
orthorhombic multiferroics, namely GdMnOz, TbMnOs 
0, Euo.75Yo.2^n03 DyMnOs M, l^MnaOg 

and TbMn205 [13|. However, their magnetic counter- 
part is still not clearly evidenced. Morevover, as optical 
techniques probe the zone center, the shape of their dis- 
persion and thus the underlying mechanism responsible 
for the hybridization, is still unknown. In this letter, we 
report polarized inelastic neutron scattering experiments 
performed on the particular case of hegagonal YMnOs 
to unravel the spin-lattice dual nature of these hybrid 
modes. Moreover, as neutron scattering allow a global 
survey of the reciprocal space, we report the dispersion of 
these hybrid modes throughout the Brillouin zone. This 
result is discused in the framework of the dynamical mag- 
netoelectric coupling theory, where the Dzyaloshinskii- 
Moriya interaction plays a central role. 

YAInOs becomes ferroelectric below 900K, with an 
electric polarization along the c-axis, due to alternatively 
long and short yttrium-oxygen bonds (parallel to the c- 
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FIG. 1: (Color online) Dynamical magnetic and nuclear 
changes from the paramagnetic phase to the Nel ordered state. 
Figures (a) to (d) report energy scans in the SF (red la- 
belled) and NSF (blue labelled) channels at the wave vector 
Q=(0,0,6) measured in the paramagnetic phase (left column) 
and in the Neel state (right column). Red and blue lines are 
fits to the data including the spin waves and the nuclear mode 
for low temperature scans. 



axis) [13, ■ Despite a strong geometric frustration, the 
Mn spin order below the Neel temperature T/v = 75K in 
a classical triangular arrangement [l^, 17, [l3| (Curie- 
Weiss temperature 9 500K). The magnetic and electric 
order parameters strongly interact, as recently shown in 
reference [l9j , claiming the existence of a giant magneto- 
elastic coupling. 

Through the comprehensive investigation of hybrid 
spin and lattice excitations, inelastic neutron scattering 
on a triple axis spectrometer combined with longitudi- 
nal polarization analysis (LPA) [2C' 



2l| offer an efficient 



way to determine hybrid modes. Indeed, this technique 
allows to measure separately the spin-spin and nuclear- 
nuclear correlation functions in different channels (called 
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spin flip SF and non spin flip NSF respectively. It is 
worth noting that the magnetic cross section probes spin 
fluctuations perpendicular to the wave vector Q only. 
On the contrary, the nuclear cross section probes the 
atomic displacements parallel to Q. Measurements were 
carried out on the IN22 triple axis spectrometer at the 
Institute Laue Langevin (ILL, G renoble, France). The 
sample was aligned in the scattering plane (100), (001) 
such that momentum t ransfer of the form Q=(II,0,L) 
in units of reciprocal lattice wave vectors were accessible 
and mounted into a n ILL- type orange cryostat (1.5K- 
300K). All data were obtained with a fixed final wave 
vector of 2.662A~^ providing an energy resolution less 
than ImeV. Heussler crystals were used as analyzer and 
monochromator, together with a flipper of Mezei to re- 
verse the spin of the scattered neutrons. The elastic and 
inelastic measurement of the polarisation efficiency (flip- 
ping ratio), as determined from different magnetic Bragg 
peaks (100), (105), (003) and from the magnon signal at 
Q=(-0.4 6) and (1 1), was of about 16 (elastic) and 
14 (inelastic signal). The amplitude of the expected leak 
from the SF channel to the NSF channel is therefore 
less than 7% of the SF intensity. Using an unpolarized 
beam (PG monochromator), we measured a polarization 
parallel and perpendicular in- and out-of-plane to the 
wave vector less than 10~^ excluding the presence of chi- 
ral terms and of the nuclear-magnetic interference terms 

Energy scans at the wave vector Q=(0,0,6) in both 
channels measured above and below Tjv are depicted on 
the panel B of the Figure [TJ On the one hand, the strong 
magnetic quasi-elastic signal observed at 150K in the SF 
channel confirms the presence of strong spin-spin corre- 
lations in the paramagnetic phase arising from the geo- 
metrically frustrated Mn moments 3i 13 • Below T/v, 3 
almost doubly degenerate spin wave modes are known 
to rise up [l^, 122, l23|. At the zone centre, they are 



characterized by a spin gap Asg, typical of a magnetic 
anisotropy. Upon cooling below T^r, the magnetic long 
range order develops and Asg gradually shifts to higher 
values (Figure 2c), reaching its maximum of about 5.3 
meV around 40K. On the other hand, the NSF inten- 
sity shows at 150 K a quasi-elastic signal as well. At 
1.5K, the NSF data demonstrate the emergence of an 
additional inelastic nuclear mode. The energy of this 
mode coincides with the spin gap Age, pointing out its 
close connection with the spin subsystem. To understand 
the temperature dependence of both quasi-elastic and in- 
elastic modes, NSF energy scans have been performed 
at different temperatures around (Figure [5^). From 
these NSF data, we extract the lattice susceptibilities, 
presented on Figure [21d. This analysis shows that the ad- 
ditional mode rises upon cooling on the top of the quasi- 
elastic signal as a supplementary intensity. Below 40K, 
the energy of this nuclear mode is found to follow the spin 
gap Asg (Figure Ht). Above 40K, its presence becomes 



however hard to distinguish. To overcome this difficulty, 
we measured the NSF intensity at Q= (0,0,6) for sampling 
temperatures at 2.2 meV, well below the low tempera- 
ture energy of the nuclear mode (Figure [2ji) . According 
to the above analysis, the quasi-elastic signal is expected 
to give a contribution proportional to the standard de- 
tailed balance factor (grey line). Figure [2] shows however 
subtle deviations from it in the intermediate region sur- 
rounding Tjv (blue line) , that we attribute to the raising 
of the nuclear mode starting from the lowest energies. 
These measurements allow estimating the energy of the 
nuclear mode for temperatures higher than 40K. Finally, 
we conclude that both spin gap a nd nuclear mode en- 
ergies follow the same temperature dependence from 1.5 
K till T/v. Since the nuclear cross section probes fluctua- 
tions along Q=(0,0,6), which is parallel to the c-axis, we 
attribute the nuclear quasi-elastic signal to relaxational 
vibrations along this particular direction. The nuclear 
mode corresponds to collective vibrations along the same 
c-axis. Moreover, since it is found at a zone center, it cor- 
responds to vibrations within the decoration of the unit 
cell. Now the question is: how do such internal motions 
propagate through the crystal? 

To address this question, it is essential to determine 
the dispersion of the additional nuclear mode, by repeat- 
ing the same measurements for different Q=(H,0,6) with 
varying H (Figure [3|). We first determine the disper- 
sion of the modes which are trivially expected, namely 
the spin wave modes (red points in SF channel on Fig- 
ure [3]) and the transverse acoustic phonon22 (denoted 
by P on Figure [3]). Amazingly, the NSF data demon- 
strate unambiguously the existence of additional nuclear 
modes all along the low energy spin wave dispersions from 
Q=(0,0,6) to Q=(0.25,0,6) (denoted by HM on FigurelS]). 
Their intensity is found to decrease as the wave vector 
increases. For instance, at Q=(0.325,0,6), the nuclear 
intensity measured in the NSF channel becomes compa- 
rable to the expected leak of the intensity from the SF 
channel. Figure[3^ evidences the fact that the dispersions 
of these nuclear modes match the spin wave dispersions 
for a wide range of wave vectors. It therefore appears 
obvious that they must be attributed to collective vibra- 
tions within the decoration that propagate through the 
crystal in the Neel state by hybridizing with the spin 
waves. 

The discovery of a strong mixing down to the zone 
centre 11=0 shows that the long range properties of the 
system are affected. In explaining these findings, we thus 
propose to consider a coupling of the spin subsystem with 
atomic displacements within the unit cell [l9[. This in 
turn implies a hybridization mechanism with an optical 
phonon, as in the dynamical magnetoelectric coupling 
theory developed for orthorhombic RMnO^ 0, H, [IS]- 
In this scenario, the spin current Jij ~ Si x Sj (de- 
fined for neighbouring manganese spins Si and Sj sit- 
ting at distance r^) plays a crucial role. Owing to the 
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FIG. 2: (Color online) Nuclear mode at the zone center in 
the Nel phase, (a) Energy scans in the NSF channels at the 
wave vector Q= (0,0,6) for different temperatures above and 
below Tjv. Below Tjv, the additional nuclear excitation starts 
to form from the 40K (blue line) . (b) Background subtracted 
and detailed balance factor divided neutron intensity of a), (c) 
Temperature dependences of the spin gap (red points, mea- 
sured with unpolarized neutrons) and of the nuclear mode 
(blue points) positions, (d) Temperature dependence of the 
NSF inelastic intensity at 2.75 meV. The gray line corresponds 
to the detailed balance factor. 



spiral magnetic ordering typical of these materials, Jy 
acts, via the inverse Dzyaloshinskii-Moryia coupling, as 
a force pushing the oxygen atom located between adja- 
cent manganese off the Mn-Mn bond. This mechanism 
gives rise at T/y to an electric polarization P = rtj x 
lying within the spiral plane and perpendicular to r^. 
The multiferroic Goldstone mode is predicted to be a hy- 
brid mode, rising at AsGj and made of a mixing between 
the optical phonon associated with the oxygen displace- 
ment, and the spin wave mode involving spin fluctuations 
out of the spiral plane. At first glance, this mechanism 
would not hold in the hexagonal case, since due to the 
triangular symmetry, the resulting magnetic force expe- 
rienced by oxygen atoms is zero: indeed, each oxygen ion 
is located at the centre of a triangle formed by 3 neigh- 
bouring Mn ions. The system can however benefit from 



FIG. 3: (Color online) Nuclear and magnetic energy disper- 
sions. Energy scans taken at 1.5K for different wave vectors 
Q=(H,0,6) in both SF (fuU points and red fits) and NSF (open 
points and blue fits) channels. The H component is indicated 
in the top left corner of each panel. The dotted line indi- 
cates the position of the nuclear mode (see text). The lowest 
right panel depicts the dispersions of the different low energy 
excitations : measured (filled red points) and calculated (dot- 
ted red lines) spin waves, transverse acoustic phonon6 (blue 
square) and hybrid modes (open blue circles). 



the Dzyaloshinskii-Moryia interaction by spontaneously 
moving the oxygen atoms along the c-axis. This distor- 
tion is expected to create an electric polarization parallel 
to Jij and is accompanied by a slight rotation of the 
spins towards the same direction. In that case, the cou- 
pled spin and atomic motions look like those of the ribs 
of an umbrella that would be put up or down (Figure [?]). 
In close analogy with the orthorhombic case, hybridized 
spin-lattice Goldstone modes are expected at A513, in 
agreement with the present results. 

To test the validity of this scenario, several predictions 
have to be further examined. First, the magnetic struc- 
ture should be characterized by a tiny ferromagnetic mo- 
ment superimposed on the AF structure, as sketched in 
the right panel of Figure 01 In that case, specific Bragg 
peaks should exhibit a specific contribution at T/y- Cal- 
culations of the structure factor show that this effect is 
best observed for Bragg peaks with a forbidden AF mag- 
netic intensity and a very weak nuclear one. The (2-11) 
Bragg peak fulfils these conditions. As shown on Figure 
m its intensity increases below T^r, and this is a good 
indication for the validity of this scenario. Next, as the 
high temperature ferroelectric distortion is mainly due to 
atomic displacements parallel to the c-axis, the oxygen 
displacements proposed in this umbrella scenario should 
result at Ty in a slight change of the ferroelectric mo- 
ment. Evidence for such an evolution has been recently 
reported by Lee et al [l^, thanks to high resolution X 
rays and neutrons diffractions measurements. This result 
is another argument supporting our interpretation. 

In conclusion, polarized inelastic neutron scattering ex- 
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FIG. 4: (Color online) Umbrella scenario. Left figure (a) 
displays the temperature dependence of the (100) antifero- 
magnetic bragg peak and of the (2-11) Bragg peak (unpolar- 
ized neutron) .The umbrella mechanism discussed in the text 
is sketched on the right panel (b) . It evidences the existence of 
a small ferromagnetic component perpendicular to the Man- 
ganese (Mn)-Oxygen (O) layers when the oxygens are moving 
out of the plane. 

periments demonstrate the existence of hybrid spin and 
lattice low energy modes in YMnOs, that can be consid- 
ered as Goldstone modes of the multiferroic phase. The 
neutron polarization analysis directly shows their hybrid 
nature, revealing both spin and structural counterparts. 
The mechanism responsible for this hybridization could 
be the Dzyaloshinskii-interaction, in close analogy with 
the model recently proposed for orthorhombic multifer- 
roic materials. 
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